Hadron Physics and Lattice QCD 



Constantia Alexandrou 

Department of Physics, University of Cyprus, P.O. Box 20537, 1678 Nicosia, Cyprus and 
Computational-based Science and Technology Research Center, The Cyprus Institute, P.O. Box 

27456, CY-1645 Nicosia, Cyprus 

Abstract. A review of recent lattice QCD hadron structure calculations is presented. Important 
hadronic properties such as the axial charge and spin content of the nucleon, as well as, the mass 
and axial charge of hyperons and charmed baryons are discussed. 
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INTRODUCTION 

The recent progress in the numerical simulation of the fundamental theory of the strong 
interactions, Quantum Chromodynamics (QCD), has been remarkable. Improvement 
in algorithms coupled with increase in computational power have enabled simulations 
to be carried out at near physical parameters of the theory. This opens up exciting 
possibilities for ab Initio reliable calculation of experimentally measured quantities as 
well as for predicting quantities that are not easily accessible to experiment. During 
the last decade, results from simulations of QCD have emerged that already provide 
essential input for a wide range of strong interaction phenomena as, for example, the 
QCD phase diagram, the structure of hadrons, nuclear forces and weak decays. In this 
presentation we focus on hadron structure calculations using state-of-the art lattice QCD 
simulations [1, 2, 3]. Understanding nucleon structure from first principles is considered 
a milestone of hadronic physics and a rich experimental program has been devoted to 
its study, starting with the measurements of the electromagnetic nucleon form factors 
initiated more than 50 years ago. Reproducing these key observables within the lattice 
QCD formulation is a prerequisite to obtaining reliable predictions on observables that 
explore Physics beyond the standard model. 

The starting point of lattice QCD (LQCD) is a definition of the theory on a four- 
dimensional Euclidean space-time lattice with gauge fields defined as links between ad- 
jacent lattice sites and quarks defined at each lattice site as anticommuting Grassmann 
variables belonging to the fundamental representation of SU(3) [4]. It can be simulated 
on the computer using methods analogous to those used in Statistical Mechanics allow- 
ing calculation of matrix elements of any operator between hadronic states in terms of 
the fundamental quark and gluon degrees of freedom. The discretization the QCD La- 
grangian is not unique and this is reflected in the variety of lattice actions of QCD that 
are in use today. It is customary to classify them according to the discretization of the 
fermionic part of the Lagrangian, the most common being improved Wilson-type such 
as Clover and twisted mass (TM) fermions, staggered fermions, domain wall (DW) and 



overlap fermions. The latter two preserve chiral symmetry but are more expensive to 
simulate. Although LQCD provides an ab Initio calculation of hadronic properties, the 
discretization of space-time and the numerical simulation on a finite volume introduce 
artifacts that may lead to systematic errors, which must be carefully investigated before 
comparing to experimental results. Since simulations become increasingly harder as the 
quark mass is decreased to its physical value, LQCD simulations were up to very re- 
cently being carried out for quark masses larger than physical. Currently, simulations 
close or even at the physical value of the pion mass are being produced and analyzed. 

RECENT RESULTS 
Hadron masses 

Masses of low-lying hadrons are considered benchmark quantities for LQCD. They 
are computed by evaluating the vacuum expectation value of two-point functions. Recent 
results using Clover [5] and TM fermion actions [6] with a careful analysis of lattice 
systematics have been obtained showing agreement with the experimental values. The 
QCDSF-UKQCD [7] and PACS-CS collaborations [8] have also performed spectrum 
studies using Clover fermions. Recent results using staggered [9], domain wall fermions 
(DWF) [10], and a mixed action approach with staggered sea quarks and DW valence 
quarks [11] show overall consistency [12]. This is a significant validation of LQCD. 

There has been remarkable progress recently in computing masses of excited states 
of hadrons, notably by the Hadron Spectrum Collaboration using anisotropic lattices 
with a fine lattice spacing in the temporal direction that leads to an improved resolution 
of excited states [13]. Experimental searches of charmed hadrons have received signif- 
icant attention, mainly due to the experimental observation of candidates of the doubly 
charmed baryons S+ (3520) and S++(3460) by the SELEX collaboration [14, 15, 16]. 
No evidence was found for these states by the BABAR experiment [17] and FOCUS 
Collaboration [18]. The BELLE Collaboration [19] finds S-states lower in mass, that 
can be candidates of excited states of S c but no doubly charmed cascade. Therefore, 
it is interesting to compute these masses in LQCD and compare with the experimental 
values. A mixed action approach is in general adopted with the charm valence quark 
being introduced on dynamical gauge configurations, which are produced with either 
staggered sea [20, 21, 22] or TM [23] fermions. A comparison of recent LQCD results 
on the masses of charmed baryons is provided in Fig. 1 . There is an overall agreement 
among LQCD results apart from the mass of S cc . Although for E cc a value consistent 
with the result of the SELEX experiment is obtained with TM fermions (TMF), further 
study is required to understand the results among different lattice computations in order 
to reach a final conclusion. In Fig. 1 we also compare results for the spin 3/2 charmed 
baryons [23] . There is good agreement among lattice results and with the known experi- 
mental values for E*, E* and £1*. Thus the lattice results can be taken as a prediction for 
the masses of the charmed spin-3/2 baryons S* c , Cl* c and Cl C cc- 
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FIGURE 1. Masses of charmed baryons with spin 1/2 (left) and spin 3/2 (right) computed within 
lattice QCD. The experimental values are shown by the horizontal lines. With red filled circles, we show 
TMF results extrapolated to the physical pion mass whereas open circles shows the results obtained at 
m K = 260 MeV [23]. Results obtained using a number of hybrid actions with staggered sea quarks are 
from Refs. [21] (open blue squares), [20] (open green triangles) and [22] (filled magenta triangles). 



Hadron From factors 

Calculation of hadron matrix elements of the form (hf(pf)\0\hi(pi)) requires both 
the evaluation of two- and three-point functions [2]. 

Nucleon axial charge: There is a number of LQCD calculations of the nucleon axial 
charge gA, which is considered a benchmark quantity for form factor calculations within 
LQCD. The reason for this is that it is both well known experimentally and it can 
be determined at momentum transfer squared q 2 = (pf — pt) 2 = with no ambiguity 
associated with having to fit the g 2 -dependence of the form factor (FF), as for example, 
in the case of the anomalous magnetic moment where one needs to fit the small in- 
dependence of the magnetic FF In addition, only the connected diagram where the axial 
current couples to a valence quark as shown in the diagram, 
contributes. The computational cost is about twice that of 
a two-point function involved in the calculation of hadron ■''> 
masses. 

In Fig. 2 we show recent LQCD results using TMF, Clover fermions, DWF and a hy- 
brid action of DWF on a staggered sea, all of which are renormalized non-perturbatively. 
As can be seen, there is a nice agreement among different lattice discretizations and no 
significant dependence on the quark mass down to about m n = 250 MeV. Clover results 
with N f = 2 + 1 at m n ~ 150 MeV [24] and N/ = 2at m K ~ 1 80 MeV [25] underestimate 
the value of gA- The origin of this discrepancy is being investigated. 
A possible explanation for the observed discrepancy could come from excited states 
contamination. A dedicated high accuracy study of excited state contributions showed 
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FIGURE 2. Left: Recent LQCD results on g A using Nf = 2 TMF [26] (filled circles), Nf = 2 + 1 + 1 
TMF (filled squares), Nf = 2+1 DWF [27] (crosses), Nf = 2 + 1 hybrid action [28] (filled triangles) 
Nf = 2 Clover [25] (square with cross) and Nf = 2+1 Clover [24] (open circles). The experimental value 
is shown by the asterisk. Right: Study of excited state contributions on gA for various time source-sink 
separations t s . The light blue band is the result of summation method. The upper plateau is for the isovector 
axial current whereas the lower for the isoscalar. 



that for m % ~ 400 MeV there is no significant effect on the value extracted for gA [29] . 
This is illustrate in Fig. 2 where the time separation between sink and source is varied 
from t s = 4a to t s = 12a. As can been seen, the same plateau value is obtained for all 
time separations. The same value is also obtained from another approach that involves 
summing the time where the axial current couples to the quark inside the nucleon. This 
demonstrates that excited state contamination is negligible at least for m K ~ 400 MeV. 
An assessment of volume and cut-off effects was also carried out [26] indicating that, 
for pion masses larger than m n ~ 300 MeV, volume and discretization errors are small 
compared to the uncertainty in the chiral extrapolation. This uncertainty is, however, 
eliminated by the recent result at almost physical pion mass, which, if confirmed by 
further studies, confronts us with a challenge. 

Axial charge for hyperons and charmed baryons: The axial couplings of hyperons 
are either less well measured or not known experimentally. One relies on theoretical 
estimates, which can have large uncertainties. These couplings are phenomenologically 
important parameters within effective field theory descriptions. Using LQCD one can 
evaluate these axial couplings using the techniques employed in the case of the nucleon 
gA, i.e. they are determined by the connected part of the appropriate hadron matrix 
element of the axial vector current at zero momentum transfer. 

A first lattice QCD calculation of the axial charge of spin- 1/2 hyperons was carried 
out using staggered sea quarks and DW valence quarks [30]. Using Nf = 2 + 1 + 1 
TMF the axial charges of hyperons and charmed baryons are calculated using the 
fixed current method that enables the computations of the axial charge of all hadrons 
with one sequential inversion, i.e. with the same number of inversions as for the case 
of the nucleon axial charge. The SU(3) symmetry relation g^ = F + Z), = 2F, 
g^ = —D + F leading to g% — g^ + g^ =0 is shown as a function of the breaking 
parameter x = (m 2 K — mfyjiAn f%) in Fig. 3 for both TMF and hybrid results. The 
SU(3) symmetry breaking <5 S tj(3) = g 1 ^ — g\ + nas an * 2 -dependence yielding at the 
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FIGURE 3. Left: <5 SU(3) = g% - g\ + g% versus x for TMF (red filled circles) and hybrid (blue filled 
squares) results [30]. The line is a fit to all data whereas the blue band to the TMF results. Right: The 
axial charge of Q. and Q.® versus m^. 



physical point <5 SU (3) ~ 15%. Fist results for the decuplet axial charges g%, g^ and gjf 
indicate an even smaller SU(3) symmetry breaking. In Fig. 3 we show first results using 
Nf — 2+1 + 1 TMF for the Q. axial charge g^ as a function of m%. 
Nucleon Dirac and Pauli isovector radii and anomalous magnetic moment: The nu- 

cleon Dirac and Pauli radii can be determined from rf 2 = — j^jj^ ~^\q 2 =o wnere me 
form factors F\ and are extracted from nucleon matrix element of the electromagnetic 
current: (N(p f )\f(0)\N( Pl ))=u N (p f ) [j^F^q 2 ) + l °^F 2 (q 2 )] u N ( Pi ).ltis custom- 
ary to use a dipole Ansatz to fit the g 2 -dependence of F\ and F%. The anomalous magnetic 

phys 

moment is then given by ^(0)— \^ in Bohr magnetons. A number of LQCD collabo- 

m N 

rations have computed the isovector nucleon Dirac and Pauli form factors F\ and F 2 , 
which in the isospin limit, gives the FFs of the proton minus that of the neutron. These 
isovector quantities only receive connected contributions and are thus computationally 
straight forward to calculate. 
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FIGURE 4. Dirac (left), Pauli (middle) rms radii and the nucleon anomalous magnetic moment (right) 
using Nf = 2 [31] and Nf = 2+1 + 1 [32] TMF, Nf = 2 Clover fermions [33], Nf = 2+1 DWF [34, 27] 
and a hybrid action of staggered sea and DW valence quarks [28]. 

In Fig. 4 we collect the most recent LQCD results on the isovector Dirac and Pauli 
radii, and anomalous magnetic moment. As can be seen, LQCD results are in agreement 
even before taking the continuum limit, which indicates that cut-off effects are small. 
Recent results at m n ~ 180 MeV [33] show an increase towards the physical value. The 
less rapid fall-off of baryon FFs with q 2 , responsible for the smaller values of the radii, 



is a common feature in current lattice QCD calculations that is under investigation. 



Nucleon spin 

In order to extract information on the spin content of the nucleon one needs to 
evaluate the isoscalar moments A 1 ^ and B 1 2 'q d since the total spin of a quark in the 
nucleon is given by ]i = I (A q 2Q +B q 2Q ) . The total spin can be further decomposed into 
its orbital angular momentum L q and its spin component AE 9 as J q = ^AT, q + L q . The 

spin carried by the u- and d- quarks is determined using AT, u+d = A"q = . In order to 
evaluate the isoscalar quantities one would need the disconnected contributions. These 
are notoriously difficult to calculate and they are neglected in most current evaluations. 
A recent calculation of the quark spin, which included the disconnected contributions 
was carried out using Clover fermions giving Au + Ad + As = 0.45 (4)(9) with As = 
-0.020(10) (4) at jU = \flA GeV [35]. This calculation shows that the strange quark 
contribution to the nucleon spin is small, a result consistent with COMPASS [36] and 
HERMES data [37]. 
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FIGURE 5. Left: (x) u - d using N f = 2 [38] and Nf = 2 + 1 + 1 [32] TMF, N f = 2 clover fermions [25, 
39], hybrid [28] and DWF [40]. The physical point shown by the asterisk is from Ref. [41]. Right: Quark 
contribution to the nucleon spin and angular momentum using Nf = 2 and Nf = 2+1 + 1 TMF. 

Neglecting the disconnected contributions to the nucleon spin, one can extract J q and 
AL q using current lattice QCD data. In Fig. 5 we compare recent LQCD results on the 
isovector momentum fraction A20 = {x) u -d m me MS scheme at a scale /! = 2 GeV. 
There is an overall agreement among lattice QCD results and therefore we use TMF 
results to study the spin content. In Fig. 5 we show the pion mass dependence of the 
isoscalar spin AE and angular momentum L carried by the u- and d-quarks and in Fig. 6 
we show J u and J d as well as the angular momentum and spin for the u and d quarks. 
The expressions for the isovector and isoscalar J and AE used in the chiral extrapolations 
are 

J- d = <m 2 , + ^(l- ? -^((2gi + l)log"|+2gi) 



ju+d = 45 + afm 2 +fl /5 (1 _ 3g 2_^_ log ^| ) 

2 2 

AE = a x m K + a (l-3g A 2 log^-) (1) 

At the physical pion mass we obtain a total spin of J u ~ 1/4 and 7 rf ~ for the u- and 
d-quarks, respectively. These results are in agreement with those by the LHPC [28]. 
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FIGURE 6. Spin and angular momentum for u- and d- quarks, for Nf = 2 [38] (filled red, green and 
magenta symbols) and Nf = 2 + 1 + 1 [32] (filled blue square and black cross) TMF. 



CONCLUSIONS 

Recent simulations of lattice QCD successfully reproduce the low-lying baryon spec- 
trum using different discretization schemes paving the way for studying nucleon struc- 
ture. Key observables such as the nucleon axial charge, rms radii and anomalous mag- 
netic moment, spin and angular momenta carried by quarks are being investigated using 
simulations with parameters close to their physical values. Understanding well measured 
quantities, such as g A , is crucial in order to reliably predict other less known observables. 
Simulations currently produced with pion masses below 200 MeV, combined with a de- 
tailed study of lattice systematics, are expected to shed light both on the origin of the 
observed discrepancies and make reliable predictions on other less well known observ- 
ables. 
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